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Using harmonic emission from circularly polarized undulator is a procedure that is normally em- 
ployed on synchrotron beamlines in order to extend the covered spectral range. A similar capability 
is likewise beneficial for next generation free-electron lasers. In this paper, we perform a first quanti- 
tative experimental analysis of the angular distribution of free-electron laser harmonic emission from 
helical undulators. Experimental results are compared to those obtained by means of a theoretical 
model based on the paraxial solution of Maxwell's equations. 

PACS numbers: 



I. INTRODUCTION 

Variable polarization is one of the most attractive fea- 
tures of light sources based on electron beams propagat- 
ing through APPLE-type undulators. 

Such a feature is particularly attractive for those ap- 
plications aiming to investigate the local symmetry of 
a given system, especially when the symmetry can be 
ascribed to the lattice geometry (as in the case of crys- 
tals), to the chirality of a molecule or to a possible lo- 
cal atomic magnetic moment. In fact, several spectro- 
scopic methods rely on the opportunity to chose a well 
defined polarization state of the light. Spin [l[ and angu- 
lar resolved photo-electron spectroscopy, x-ray linear and 
circular dichroism or resonant scattering of polarized x- 
rays @, H| are only few examples. For this reason, light 
polarization-dependent spectroscopy using standard syn- 
chrotron radiation has become a very powerful tool to 
study the electronic and magnetic properties of matter 

aff- 

Thanks to a dramatic increase of photon peak- 
brilliance with respect to conventional sources, as well 
as to the possibility of controlling both the temporal du- 
ration and the spectral bandwidth of the produced light 
pulses, free-electron lasers (FELs) will allow completely 
new studies in the above mentioned research fields [fj. 
Such studies will also take advantage of well-defined and 
easily tunable light polarization. For the preparation 
of the scientific case of single-pass FELs, a fundamen- 
tal question concerns the possibility of generating signifi- 
cant photon flux at higher harmonics of the fundamental 
wavelength, while maintaining polarization ductility. 

In a FEL, harmonic generation is driven by the 
electron-beam interaction with a radiation at a given 
fundamental wavelength, in the presence of the static 
and periodic magnetic field generated by an undulator 
[7h14|. Such interaction produces a significant modula- 
tion (called bunching) of the electron-beam density at the 



fundamental wavelength, and its harmonics. Harmonic 
bunching is the (nonlinear) source of FEL harmonic emis- 
sion. The latter is orders of magnitudes stronger than the 
spontaneous incoherent emission and provides the pos- 
sibility to extend the FEL tuning range towards short 
wavelengths. 

In the simplest possible configuration, FEL harmonics 
are produced in the same undulator where radiation at 
the fundamental wavelength is generated. This method 
is the standard one when the FEL light is obtained from 
Self Amplified Spontaneous Emission (SASE) [l5|, and 
is normally referred to as nonlinear harmonic generation 
(NHG). 

An alternative approach to NHG for producing coher- 
ent harmonic emission relies upon the bunching created 
in an undulator, called the modulator, for generating co- 
herent emission in a subsequent undulator, called the ra- 
diator, tuned at one of the harmonics of the modulator 
fTfjl This technique is usually called coherent har- 
monic generation (CHG). The implementation of CHG 
requires the use of an external coherent source, acting as 
a "seed" , which is normally provided by one of the har- 
monics of a Ti: sapphire laser. These harmonics can be 
produced either using a solid-state crystal, or a gas jet 
[lq . The seed-electron interaction in the modulator is 
necessary to initiate the bunching process. Bunching is 
produced at all (odd and even) harmonics. As a conse- 
quence, in the radiator, coherent emission may occur at 
any selected harmonics, no matter the radiator polariza- 
tion. The situation changes when NHG is considered. It 
is well known that, in case of planar undulators, on-axis 
NHG occurs only at odd harmonics, while even harmon- 
ics are present only off-axis 0. In [l9| it has been ex- 
perimentally demonstrated that, as in the case of sponta- 
neous emission, the NHG signal generated by helical un- 
dulators is distributed off-axis. A quantitative estimate 
of the off-axis coherent harmonic flux is, however, still 
missing. This lack of information is due to the absence, 
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FIG. 1: Reference system and symbols used to develop the 
theoretical model. 



so far, of an experimental characterization of the angular 
distribution of harmonic emission and, at the same time, 
to the difficulty of modelling the coherent emission from 
a modulated electron beam in the far zone. Recently, 
an analytical model has been proposed [20U23I ] . relying 
on the paraxial solution of Maxwell's equations. Such 
a model confirms what has been found in i.e. no 
on-axis NHG in helical undulators, and provide a quan- 
titative prediction of the far-field angular harmonic dis- 
tribution. 

In this paper, we perform a first quantitative experi- 
mental analysis of the angular distribution of FEL har- 
monic radiation generated by a helical undulator. In par- 
ticular, we characterize the off-axis emission in NHG con- 
figuration. Obtained results are used to benchmark the 
theoretical model proposed in poM^ . 

The paper is structured as follows. In Section [TT] we 
briefly review the theoretical model and discuss the an- 
gular characteristics of the predicted far-field harmonic 
distribution. In Section IIIII we present the experimen- 
tal device we used to perform the measurements, i.e. the 
FEL installed at the Elettra storage-ring [25|, [26[ . In Sec- 
tion IIV[ experimental data are presented and compared 
to theoretical predictions. Finally, in Section [V] we draw 
conclusions and provide perspectives for future work. 



II. THEORETICAL FRAMEWORK 

The coherent emission in the far zone from modulated 
electron beams in planar or helical undulators can be 
studied by solving the paraxial Maxwell's equation for 
the electric field in the space- frequency domain poU^ . 

Let us call the transverse electric field in the space- 
frequency domain E±(z,r±,Lu), where, as shown in Fig. 
[TJ r ']_ = xe x + ye y identifies a point on a transverse plane 
at longitudinal position z, e x and e y being unit vectors 
in the transverse x and y directions. 

The paraxial approximation is well suited for mod- 
ellingthe emission from ultra-relativistic electrons (see 
[2H HH). It relics on the assumption that the elec- 



tric field envelope E± = E± exp [— iwz/c] does not vary 
much along z on the scale of the reduced wavelength 
A = \/{2n). The wave equation for the slowly- varying 
envelope of the electric field is found to be 
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(where c is the speed of light in vacuum). The right 
hand side is determined by the trajectory of the source 
electrons, and is written in terms of the Fourier transform 

of the transverse current density, jj_(z, f±, uj), and of the 
charge density, p(z,r±,uj), which will be considered as 
given macroscopic quantities. The Fourier transform of 
the charge density can be written as 



p = -p(z,f± - r' j_(z),u>)exp 



s (z) 



(2) 



where the minus sign on the right hand side is intro- 
duced for notational convenience. In the latter expres- 
sion, the quantities r' ±(z), s (z) and v are, respec- 
tively, the transverse position, the curvilinear abscissa 
and speed of a reference electron with Lorentz factor 7 
that is injected on axis with no deflection and is guided 
by the undulator field. Such electron follows a trajectory 
r' ±(z) — r' ox e x + r' oy e v , which assumes different forms, 
namely: 



r' ±(z) = —r- [(cos(k w z) - l)e x + sin(k w z)ey\ (3) 
if the electron propagates through a helical undulator and 



r' ±(z) 



K 



(cos(k w z) - l)e x 



(4) 



if the undulator is planar. Here K = \ w eH w / (2irm e c 2 ) 
is the undulator parameter, X w = 2ir/k w being the un- 
dulator period, e the electron charge, H w the maximal 
modulus of the on-axis magnetic field of the undulator, 
and m e the rest mass of the electron. The corresponding 
velocity is given by v ±{z) = v ox e x + v oy e y . 

Since we want to discuss an FEL process, p in Eq. ([2]) 
is a slowly varying function of z on the wavelength scale 
and is peaked around the harmonics of the fundamental 
(i.e. at frequencies hu) r = 2hk w cj 2 , j z being the average 
longitudinal Lorentz factor, and h the harmonic number) 
with bandwidth Au>/(htd r ) <S f . 

Moreover, for each particle in the beam, the relative 
deviation of the particles energy from jm e c 2 is small, 
i.e. (7 — 7")/7 <S 1. Therefore, in first approximation, we 
can neglect the difference between the average transverse 
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velocity of electrons (v±) and v ±, so that j± = (v±)p — 
v ±p- 

Furthermore, we will be interested in the total power 
emitted and in the directivity diagram of the radiation 
in the far zone. We therefore introduce the far zone ap- 
proximation calling the observation angle 9 = f± /z (see 
FigfTJ), z being the distance from the middle of the un- 
dulator to the observer. We also set 9 = \9\ and take the 
limit for z D 3> L w , where L w — N W X W is the undulator 
length. 

Finally, we introduce a coherent deflection angle, ff c >, 
to describe the transverse deflection of the electron beam 
as a whole [3(| • We therefore perform the following sub- 
stitution: 



r' (z) — > &\z,ff c) ) = ? ' {z) + rf c) z 



(5) 



It follows that also v (z) should be substituted with 
v(z,fft°)), and that ^ z {z,ff c ^) is now a function of both 
z and ff c '. With these prescriptions, we obtain the fol- 
lowing solution of Eq. JT]) : 



its h-th integer multiple. We specify how close the fre- 
quency u> is to the h-th harmonic by defining a detuning 
parameter Ch- 



C h 



where lo = hui r 
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(10) 



and N w 3> 1, major simplifications arise because fast 
oscillating terms in powers of exp[ik w z'], which expli citly 
appear in Eq. (JB|) after substituting the trajectory [241 ] . 
effectively average to zero after integration in dz'. To 
complete the resonance approximation, we further select 
frequencies such that 
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This condition on frequencies automatically selects the 
observation angles of interest: h(9 — ff c *>) 2 ^ Un- 
der the constraint imposed by (fTT|) , independently of the 
value of K, we obtain from Eq. iR ?]) analytical expressions 
for the angular spectral fluxes [32|. These expressions de- 
scribe the coherent emission in the far zone from modu- 
lated electron beams in planar or helical undulators, and 
strongly depend on the Fresnel parameter N XtV defined 
as 



Substitution of the equations of motion Eq. ((3]) or Eq. 
dU, and of an appropriate expression for p into Eq. ((6]), 
yields the electric field in the far zone. We will consider 
a simple model where 
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<j XtV being the electron beam rms transverse size, and 
/(w) the Fourier transform of the bunch temporal profile: 
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Here bh is the bunching at the h-th harmonic [3l| . 

This model describes a bunched electron beam with 
Gaussian transverse profile, with current I(t), longitudi- 
nally modulated at frequency hui r . We are interested in 
a frequency range near the fundamental harmonic, uj r , or 



cL w 



(12) 



In the helical case, the spectral flux at the h-th harmonic 
is given by: 
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Figures [5] and [3] show the angular spectral flux as ob- 
tained from Eq. (Tl3| for the CHG and NHG, respectively. 
For the calculation, use has been made of the parameters 
reported in Table U corresponding to the experimental 
setup described in the following Section. 



FIG. 2: Angular spectral flux in CHG mode, as predicted 
by Eq. (|13l) . For the calculation, use has been made of the 
parameters reported in Tabled] corresponding to the case of 
the Elettra storage-ring FEL. In this case, light is emitted 
at the radiator fundamental wavelength (i.e., h = 1, corre- 
sponding to 260 nm), which is the third harmonic of the seed 
wavelength (780 nm). . 



FIG. 3: Angular spectral flux in NHG mode, as predicted by 
Eq. (|13|l . For the calculation, use has been made of param- 
eters reported in Table U corresponding to the case of the 
Elettra storage-ring FEL. In this case, light is emitted at the 
second harmonic of the radiator fundamental wavelength (i.e., 
h = 2, corresponding to 260 nm), which is the third harmonic 
of the seed wavelength (780 nm). 
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III. EXPERIMENTAL SETUP 

The predictions of the model presented in the previ- 
ous Section have been tested on the Elettra storage-ring 
FEL, whose scheme is shown in Fig. 01 The setup, which 
is based on two independent APPLE II undulators sep- 
arated by a magnetic chicane, is very flexible: depend- 
ing on the undulator configuration, both NHG and CHG 
schemes can be implemented. The laser pulse emitted 
by a Tksapphire laser (fundamental wavelength at 780 
nm) interacts with the electron bunch within the mod- 
ulator, inducing an energy modulation in the electron 
beam. After the conversion of the energy modulation 
into spatial bunching, occurring into the magnetic chi- 
cane, the bunch enters the radiator and emits coherent 
radiation at one of the harmonics of the seed laser (the 
third one, i.e. 260 nm, for the experiments presented in 
this paper). The harmonic radiation is first monochro- 
matized and then acquired on a digital oscilloscope by 
means of a photomultiplier tube (PMT) , see Fig. 01 The 
latter is fast enough to resolve the pulsed dynamics of the 
electron bunch photoemission (nanonosecond scale), but 
does not allow to resolve the photon pulse shape (tens 
of ps). A diaphragm of 1.4 mm is placed downstream 
the radiator, before the detector. The diaphragm is used 
to select only the portion of radiation emitted within a 
small cone close to the undulator axis (angular accep- 
tance of about 100 /irad). By moving the diaphragm, it 
is possible to characterize the angular distribution of the 
produced coherent radiation. 
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FIG. 4: Schematic layout of the Elettra storage-ring FEL. 



For the reported experiments, Elettra was operated in 
single-bunch mode at 0.75 GeV (7 ~ 1468), the peak 
bunch current and its revolution period being, respec- 
tively, ~ 10 A and 864 ns. The seed repetition rate is 
300 Hz. The laser pulse and the electron bunch duration 
are, respectively, 120 fs (rms) and about 35 ps (rms). A 
stable laser-electron synchronization with a jitter of few 
ps guarantees an efficient and stable seeding process, 
with shot-to-shot fluctuations of the order of few percent. 
The main experimental parameters, also used to obtain 
the theoretical angular spectral flux reported in Figs. [5] 
andfSl are reported in TablcHJ More detailed information 
about the experimental setup, as well as about the per- 
formance of the Elettra storage-ring FEL, can be found 
in El. 



IV. RESULTS 

In this Section we present the experimental data and 
perform a quantitative comparison between experiments 



5 



TABLE I: Main experimental parameters. 



Electron beam 



Peak current 


10 


A 


Normaized energy (7) 


1468 




Bunch length (rms) 


35 


ps 


Relative energy spread 


~ 0.1 


% 


Bunching (bh) 


0.6 




Horizontal transverse size (rms) (a x ) 


200 


jim 


Vertical transverse size (rms) (a y ) 


10 


;im 


Seed laser 


Power 


~ 10 


GW 


Wavelength 


780 


nm 


Pulse length (rms) 


~ 120 


fs 


Polarization 


Horiz., Circ. 




Radiator 


Periodicity 


0.10 


m 


Number of periods (N w ) 

2irc 

Resonant wavelength ( ), CHG mode 


20 




260 


11111 



Undulator parameter (K) in CHG mode 3.19 

27TC 

Resonant wavelength ( ), NHG mode 520 nm 

Undulator parameter (K) in NHG mode 4.63 

Polarization Horiz, Circ. 




-400 -200 200 400 

6 y (nrad) 



FIG. 5: Intensity of the CHG signal at 260 nm, third har- 
monic of the seed wavelength, as a function of the angle 8 y 
(vertical direction, see Fig. [2} between the undulator axis and 
the line from the center of the ondulator to the position of the 
diaphragm. Open circles: measurements (right vertical scale) . 
Continuous and dotted lines: theoretical predictions (left ver- 
tical scale) obtained using Eq. (|13[) and the set of parameters 
reported in Table [T] Dotted line: integration of the curve 
shown in Fig[2] along 6 y , over an angle of 100 /irad around 
9 X = 0. In plotting experimental data, a correction factor has 
been used to take into account (slight) current variations oc- 
curred in the time period between successive measurements. 



and theory. 

A first set of measurements has been carried out in 
CHG configuration, tuning the radiator at 260 nm (third 
harmonic of the seed wavelength), in circular polariza- 
tion. The intensity detected by the PMT, as a function 
of 8 y , is represented by the open circles in Fig. [5] As 
expected (see Fig. [2]), the intensity has a maximum in 
correspondence of the undulator axis {9 y = 0). 

Due to some mechanical limitations intrinsic to our 
device, the diaphragm can be moved only downstream 
along the vertical direction. As a consequence, no mea- 
surements have been performed along positive #j,'s, nor 
along the 9 X direction. Error bars represent the re- 
sult of several independent measurements. Fluctuations 
are mainly due to FEL instabilities and to possible un- 
desired slight shifts of the electron-beam trajectory in 
the time period between successive measurements. Data 
have been acquired at quite low average electron-beam 
current (< 1 mA), that is in a regime characterized by 
relatively long beam lifetimes (order of hours). How- 
ever, since the harmonic intensity varies as the square of 
the number of charges involved into the emission process 
pfjj . even small current variations may reflect in signif- 
icant intensity fluctuations. In order to take this effect 
into account, a correction factor has been applied to ex- 
perimental data. In Fig. [5j continuous and dotted lines 
represent the theoretical predictions obtained using Eq. 
(|To| and the parameters reported in Table |TJ The dotted 
line has been obtained by integrating the curve shown 



in Fig. along 9 y , over an angle of 100 ^rad (the di- 
aphragm acceptance) around 9 X = (corresponding to a 
perfect horizontal alignment between the electro n beam 
and the undulator axis). The continuous line has been 
instead obtained by integrating the curve shown in Fig. [2] 
along y , over an angle of 100 /irad around 9 X = 80 /trad. 
The latter value corresponds to the best possible exper- 
imental accuracy in controlling the horizontal alignment 
between the electron beam and the undulator axis. 

The same (experimental and theoretical) analysis has 
been carried out in NHG configuration. In order to re- 
produce as much as possible the experimental conditions 
of the CHG configuration, the radiator has been tuned at 
520 nm (in circular polarization). This avoids any seed- 
electron interaction inside the radiator and, as a conse- 
quence, prevents any undesired degradation of the beam 
quality (i.e., an increase of the beam energy spread). 

Since 520 nm is not a harmonic of the input seed wave- 
length, no bunching is produced in the modulator at such 
wavelength. As a consequence, there is no coherent emis- 
sion at the fundamental wavelength of the radiator. How- 
ever, one obtains coherent emission through NHG at 260 
nm and 130 nm (second and fourth harmonics, respec- 
tively, of the radiator wavelength). As for the case of 
CHG, the intensity of the coherent emission at 260 nm 
is detected using the PMT. Results are reported in Fig. 
O As expected (see Fig. [3J, due to the finite aperture 
of the diaphragm, NHG emission is not exactly zero on 
axis. 

In order to check the agreement between measurements 
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FIG. 6: Intensity of NHG at 260 nm, second harmonic of the 
radiator fundamental wavelength, as a function of the angle 
By. Open circles: measurements (right vertical scale). Con- 
tinuous and dotted lines: theoretical predictions (left vertical 
scale) obtained using Eq. (| 13p and the set of parameters re- 
ported in Table HI Dotted line: cut of Fig. [3] along y , for 
9 X = 0; continuous line: integration of the curve shown in 
FigH along y , over an angle of 100 /irad around 6 X — 80 
/irad. In plotting experimental data, a correction factor has 
been used to take into account (slight) current variations oc- 
curred in the time period between successive measurements. 



and theory, one can compare the experimental and theo- 
retical ratios between the peak CHG and residual NHG 
on-axis intensities. Using the data reported in Figs. 
[5] and [5J and considering the (more realistic) theoret- 
ical continuous curves, one gets an experimental ratio 
of about 60 and a theoretical one of 140. Taking into 
account the above-mentioned experimental incertitude, 
as well as the approximations intrinsic to the theoretical 
model, the agreement can be considered satisfactory [33{ . 



cal undulators. Measurements carried out on the Elcttra 
storage-ring free-electron laser, have been compared to 
the predictions from a model based on the analytical so- 
lution of the paraxial Maxwell equations. As we have 
shown, the agreement can be considered as satisfactory. 
This result allows one to consider the theoretical model 
as a reliable tool for predicting the off-axis photon flux of 
next-generation FELs, operated in Nonlinear Harmonic 
Generation mode. 

In fact, predictions are in general quite pessimistic. 
As an example, calculations have been recently per- 
formed [13] for harmonic off-axis emission in the case 
of the FERMI@Elettra FEL [H], an FEL user facility 
presently under commissioning at Sincrotronc Trieste. In 
the case of FERMI, operating at 5 nm in circular polar- 
ization, FEL saturation is easily reached and about 3-10 12 
on-axis photons/pulse at the fundamental wavelength are 
produced. The flux reduction at the second (off-axis) har- 
monic, calculated using the theoretical model, is drastic: 
only about 10 8 photons/pulse are expected at 2.5 nm. 
An additional factor of about -10 4 is lost if considering 
the third harmonic emission at 1.33 nm [271 ]. Having 
provided an experimental validation of the theoretical 
model, we can conclude that collecting off-axis emission 
does not allow to recover a significant photon flux when 
the FEL is operated in Nonlinear Harmonic Generation 
mode, with circularly polarized undulators. 

As an alternative method to obtain significant (on- 
axis) circularly polarized harmonic radiation, one may 
consider the possibility to operate the radiator in ellip- 
tically polarized mode. Such a method gave interesting 
results in the case of spontaneous (linear) harmonic ra- 
diation (29| . As a future study, we plan to investigate its 
efficiency when applied to Nonlinear Harmonic Genera- 
tion in FEL's. 



V. CONCLUSIONS 

We performed a quantitative analysis of the angular 
distribution of Nonlinear Harmonic Generation in heli- 
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